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Abstract. Evidence for the existence of Mira Period-Luminosity rela- 
tions is reviewed including recent work on Miras with thick circumstellar 
shells. The calibration of the relation at K is discussed in detail. The na- 
ture of variables deviating from the PL relation and possible metallicity 
effects are also discussed. 
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■ 1. Introduction 

. The possibility that at least a subset of AGB variables could be used as distance 

indicators can probably be traced back to the work of Gerasimovic (1928). He 
suggested on the basis of statistical parallaxes that Mira variables showed a 
period - (visual) luminosity relation. This was confirmed and extended in later 
work (e.g. Gyllenberg 1929, 1930; Wilson & Merrill 1942; Osvald & Risley 1961; 
O-f Clayton & Feast 1969). However, the scatter in visual absolute magnitude at 

a given period, though difficult to estimate, seemed quite large (perhaps ~ 0.5 
mag). The statistical parallax results were extended to the near infrared and to 
Mb i by Robertson & Feast (1981) on the basis of extensive JHKL observations 
by Catchpole et al. (1979) of nearby Miras. 



2. The Slope of the Mira PL Relation 
2.1. The Large Magellanic Cloud 

Glass & Lloyd Evans (1981) and Glass & Feast (1982) found that a rather precise 
period-luminosity relation (a ~ 0.25mag) existed at JHK and rrn, i on the basis 
of 11 Miras in the LMC; these relations were further refined by Feast (1984) to 
include the data of Wood, Bessell & Paltoglou (1985). This early LMC work 
depended on only a few (often only one) infrared observations per star. 

An extensive programme (Glass et al. 1990) then allowed time-averaged 
JHK and bolometric magnitudes to be derived for 29 O-rich and 20 C-rich 
Miras in the LMC (Feast et al. 1989). At K the PL relation for O-Miras was 
found to be, 

Mk = — 3.47(±0.19) log P + /? (1) 

the standard deviation from this equation is only 0.13 mag. Even though some 
of this dispersion must be observational, it was found that there was some ev- 
idence for a period-luminosity-colour relation. Such a relation is of course to 
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be anticipated since one can hardly expect the Mira instability strip to have no 
width. Recently Glass & Lloyd Evans (2003) have rederived periods for some 
of these stars from MACHO data. Their results show that the above relation is 
rather stable since their rederivation for 26 stars yields a slope of —3.52 ± 0.21 
insignificantly different from the above value. In fact the small change is due to 
the fact that they leave out three stars. It could be argued that at least two of 
these might well be included as O-Miras (leading to a slope of —3.49) since their 
spectral types (Wood et al. 1985) are K? and S. Note that the designation (M), 
derived from the photometry, in Feast et al. (1989) implies only a likely O-Mira 
(i.e. K,M or S). The PL relation at K for LMC C-Miras appears to be the same 
as for the O-Miras. 

Feast et al. (1989) also found a PL(M feoZ ) relation for LMC O-Miras with a 
slope of -3.00 ± 0.24. The data of Glass & Lloyd Evans (2003) gives a slope of 
—3.06 ± 0.26 (26 stars) or —3.03 (28 stars). In the case of M^oi the slope for the 
C-Miras is distinctly different (— 1.86±0.30). However, it has never been entirely 
clear whether this difference is fundamental or due, for instance, to the method 
of deriving the bolometric magnitudes. This could introduce systematic errors 
since the spectra of the O- and C-Miras are so different (see also below). A PL 
slope in M^oi for C-Miras derived by Groenewegen & Whitelock (1996) using 
both LMC and Galactic C-Miras of known distance has a slope approximately 
midway between those for O- and C-Miras just given. 

The above discussion refers to Miras with periods shorter than about 420 
days. In the sample of Feast et al. (1989) the Miras of longer period fell above a 
linear extrapolation of the PL relations both at Mk and M^. Hughes & Wood 
(1990) chose to fit a PL relation of steeper slope for Miras of periods greater 
than 400 days. It was later suggested on the basis of approximate theoretical ar- 
guments (Willson 2000) that this apparent steeping of the PL at longer periods 
was to be expected. However, the subsets of LMC Miras so far discussed were 
all established in optical surveys. They thus are strongly biased to Miras with 
thin circumstellar shells. In our own Galaxy it is known that Miras with thick 
circumstellar shells are more frequent amongst the longer period stars. A series 
of papers (see references in Whitelock et al. 2003) described the identification of 
luminous AGB stars in the LMC from IRAS and their observation at JHKL 
and with ISO. A subset of these stars are large amplitude (i.e. Mira) variables 
and these were recently discussed by Whitelock et al. (2003). Extensive JHKL 
photometry of these stars (including data from Wood et al. 1992 and Wood 
1998) allowed the derivation of periods and, together with IRAS and ISO data, 
bolometric magnitudes. Both O-rich and C-rich objects are present in this sam- 
ple. For the C-rich objects the periods extend up to 939 days and for the O-rich 
variables up to 1393 days. Whitelock et al. discuss in detail the estimation of 
the bolometric magnitudes of these stars which can be done in various ways. 
Whilst there remains uncertainty about the precise bolometric magnitudes of 
these stars it is clear that the bulk of those with periods greater than 400 days 
fall below the steep PL relation of Hughes k. Wood (1990) and that both the 
O- and C-rich objects lie near an extrapolation of the PL(Mfe /) for O-Miras 
of shorter period. Whitelock et al. point out that the long period C-Miras in 
Magellanic Cloud clusters (Nishida et al. 2000) also fall on an extrapolation of 
this relation. Extensive mid-infrared work would be required to determine how 
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much of the scatter about the PL relation in the sample of Whitelock et al. was 
intrinsic. 

It is known that there are low amplitude semi-regular variables above the 
Mira PL relation (see e.g. Whitelock (1986), Bedding & Zijlstra (1998) Wood 
(2000)). Such stars are believed to be evolving towards the Mira PL relation. 
The nature of the long-period Miras above the PL relation has not been clear. 
However, Whitelock (see Whitelock & Feast 2000a and Whitelock et al. 2003) 
has pointed out that some at least of these Miras are known to have strong 
lithium lines in their spectra (Smith et al. 1995) and this shows them to be 
undergoing Hot-Bottom-Burning (HBB). She therefore suggests that this is the 
distinguishing characteristic of these objects. 

It will be clear from the above that considerable caution is required in using 
Miras with periods greater than 400 days as distance indicators, especially in 
optically selected samples. An example of the problems that can arise is provided 
the discovery of a 641 day, large amplitude, variable in the Local Group Irregular 
galaxy IC1613 (Kurtev et al. 2001). This lies about 0.7mag above the PL relation 
at K (for the adopted distance of this galaxy). Suspicion as to its true nature is 
aroused by the fact that it has a very early spectral type (M3) and is rather blue 
(J — K ~ 1.14) for its long period. It is possibly a HBB star or alternatively a 
relatively faint supergiant variable. 



2.2. Is the PL Slope Universal? 

Wood (1995) showed that the PL(iT) relation for SMC Miras has closely the 
same slope as that in the LMC. Miras in Galactic globular clusters fit the PL(X) 
slope from the LMC (Feast, Whitelock & Menzies 2002). In the SgrI field of the 
Galactic Bulge, Glass et al. (1995) found that the Miras (all probably O-Miras) 
fitted a PL(i^) relation of slope —3.47 ± 0.35, i.e. identical to that in the LMC. 
In this case the scatter about the relation is much higher than in the LMC (0.35 
mag compared with 0.13 mag). This difference is primarily due to the depth of 
the Bulge in the line-of-sight. 

In our own Galaxy, Glass et al. (2001) derived periods and mean K mag- 
nitudes for 406 Miras in an area 24 x 24 arcmin square around the Galactic 
centre. There is large and variable extinction in this direction. The general 
distribution of the stars in the K — log P plane (their fig 4) nevertheless suggests 
that most of these objects would fit a PL relation with the LMC slope. However, 
Oritz et al. (2002) used ISO and near IR data to derive bolometric magnitudes 
for OH/IR Miras within half a degree of the Galactic Centre. These stars with 
relatively thick circumstellar shells show a wide range in luminosity at a given 
period. Oritz et al. find that whilst many of their stars lie on a PL(M^) rela- 
tion derived from the SgrI stars of Glass et al. (1995), a considerable number lie 
below this relation. Thus there appears to be evidence, at least at long period, 
for Miras below the PL (M{, /) relation near the Galactic Centre. The matter 
cannot be considered completely settled yet in view of the considerable difficulty 
in determining the mean bolometric magnitude of the variables and also in esti- 
mating the interstellar extinction. It is worth noting that if such "subluminous" 
OH/IR Miras were present in the LMC they might well have been too faint for 
existing surveys. 
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The most striking recent work on the Mira PL slope is that of Rejkuba 
(2003) who established a well defined PL(K) relation from about 1000 Miras in 
the Giant Elliptical Galaxy Cen A. The infrared colours of these stars strongly 
suggests that they are mainly O-Miras. Most of these Miras have logP between 
2.3 and 2.7. The slope in Cen A (—3.37 ± 0.11) clearly agrees well with that in 
the LMC, -3.47 ±0.19. 

3. The Mira PL(AT) Zero-Point 

A realistic determination of the Mira PL zero-point from trigonometrical paral- 
laxes became possible following the Hipparcos astrometric mission. The Hippar- 
cos Mira parallaxes are not of sufficient accuracy to allow their individual use 
in a determination of both the slope and zero-point of a PL relation. However, 
since the results summarized above suggest that the PL(K) slope found for the 
LMC is closely the same elsewhere, we can reasonably adopt this slope as appli- 
cable to local Galactic Miras. It is then possible to use the method of reduced 
parallaxes to obtain a zero point. Whitelock & Feast (2000b) used this method 
together with extensive infrared photometry (Whitelock et al. 2000) to derive (3 
in equation 1 above for O-rich Galactic Miras. 

In Whitelock et al. (2000) (see also Whitelock 2003) it was shown that the 
shorter period Miras could be separated by their colours into two sequences, 
the "short-period-red" and "short-period-blue" stars. These two groups are 
also distinguished by their kinematics. The latter belong to the main Mira 
sequence. Whitelock et al. found evidence from the parallaxes that the "short- 
period-red" Miras were brighter at a given period than those of the main Mira 
sequence. Whilst the main Mira PL sequence is believed to represent the end 
points of AGB evolution for stars of different initial masses and/or metallicities, 
the evolutionary status of the "short-period-red" stars is still uncertain. 

Whitelock & Feast (2000b) found from Hipparcos parallaxes of O-Miras 
(omitting "short-period-red" stars) that (3 = 0.84 ± 0.14 mag. This result de- 
pends on data for 180 Miras though most of the weight resides in a relatively 
small number of stars. It has recently been suggested that the Hipparcos results 
for some Miras require a recalibration of the adopted corrections for chromatic- 
ity (Platais et al. 2003, Pourbaix et al. 2003). Knapp et al. (2003) have published 
revised parallaxes based on these new corrections. Using these new parallaxes 
Whitelock (Leiden Workshop on the Future of AGB Star Research 2003 and to 
be published) finds a zero-point of [5 = 1.04 ± 0.13 from 140 stars. The change 
in zero-point is partly due to the fact that the selection of stars is slightly dif- 
ferent since Knapp et al. do not give revised parallaxes for all the Hipparcos 
Miras. A discussion of statistical bias corrections (Feast 2002) suggests that a 
small correction of ~ 0.02mag is necessary leading to (3 = 1.06. Table 1 lists this 
value of the PL(K) zero-point together with two other Galactic determinations. 
Vlemmings et al. (2003) have obtained distances to four Mira variables using 
VLBI astrometry of OH maser spots and this leads to the value of (5 listed. The 
final figure in the table is from Miras in globular clusters with cluster distances 
based on Hipparcos parallaxes of subdwarfs (Feast et al. 2002). 

Whitelock & Feast (2000) also analysed the parallaxes of Galactic C-Miras. 
However, in view of the uncertainties, which they discuss, it would appear best 
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Table 1. The Zero-Point of the Mira PL(K) Relation 

Method (3 
Parallaxes 1.06 ± 0.13 
OH VLBI 1.01 ± 0.13 

Globular Clusters 0.93 ± 0.14 

Mean 1.00 ± 0.08 



at present to rely on the LMC result that the (unobscured) C- and O-Miras 
fit the same PL(if) relation (Feast et al. 1989). It has already been noted (see 
Whitelock et al. 2003) that within the uncertainties both the O- and C-rich Miras 
fit, the same M^-logP relation at least at long periods. 



4. Is the PL(K) Zero-Point the Same Everywhere? 

There are a number of tests that can be applied to see whether the above PL(K) 
calibration is applicable to all O-Miras. 

(1) As Table 1 shows the value of {3 from Galactic field Miras and those in 
Globular Clusters is not significantly different. 

(2) The mean value of (3 from Table 1 (1.00) when used with the results 
of Feast et al. (1989) leads to an LMC modulus of 18.48 with an uncertainty 
of about 0.1 mag or slightly less. This agrees with an LMC modulus based on 
Cepheids calibrated in a variety of manners (Feast 2003) which is 18.52 and has 
a similar small uncertainty. 

(3) The difference between the moduli of the SMC and LMC derived from 
the Mira Ph(K) relation is ~ 0.4mag (Wood 1995) which compares well with 
the difference derived from observations of the magnitudes at the red giant tip 
(0.44 mag, Cioni et al. 2000). A similar result is obtained from Miras in M^x 
(Cioni et al.2003). 

(4) The results of Glass et al. (1995) on Miras in the SgrI field of the Galactic 
Bulge together with the calibration of Table 1 yields a distance to the Galactic 
centre (Rq) of 8.5 ± 0.7kpc. This may be compared with 8.5 ± 0.5kpc derived 
from Cepheid kinematics (Feast & Whitelock 1997) and 8.0 ± 0.4kpc recently 
derived from the orbit of a star very close to the central black hole (Eisenhauser 
et al.2003). 

(5) With the Galactic calibration of Table 1, the Mira PL(K) results of 
Rejkuba (2003) for Cen A yield a distance modulus for this giant elliptical of 
28.0±0.2mag agreeing with the RGB tip distance derived by Rejkuba, 27.9 ±0.2 
mag. 

These tests indicate that there is no present evidence for any variation of 
the PL(K) zero-point at least within the observational uncertainties. 



5. Are there Metallicity Effects on the PL relation? 

It is known that at least in globular clusters there is a relation between [Fe/H] 
and Mira period (see Feast & Whitelock 2000). We then need to know whether 
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there are variations in chemical abundance at a given period between Miras in 
different systems. That such seems likely between our Galaxy and the LMC was 
shown by Glass et al. (1995). They found that the JHK colour-period relations 
were different in the SgrI field of the Galactic Bulge from those in the LMC. 
This question has been discussed by Feast (1996) and Feast & Whitelock (2000). 
The Galaxy-LMC difference can be attributed to the weakening of the infrared 
H2O band due to a relative underabundance of oxygen or an overabundance of 
carbon (locking up oxygen in CO) in the LMC. Whether this also implies an 
overall metal deficiency at a given period in LMC Miras compared with those 
in our Galaxy is not certain since the [O/Fe] ratio for these stars is not known. 
Whilst the empirical results discussed above indicate that the PL(X) relation 
holds in all the places to which it has so far been applied, the colour differences 
between the SgrI field and the LMC show that PL relations at some other colours 
may not be universal. For instance Feast & Whitelock (1999) point out that that 
there is a difference in the (K — m;> ;)-log P relations in the LMC and the SgrI 
field. This suggests that a PL(Mfe c ;) relation is more metallicity sensitive than 
a PL(K) relation, as proposed theoretically by Wood (1990). 

6. Can Other AGB Variables be used as Distance Indicators? 

Wood (2000) showed that AGB variables in the MACHO survey of the LMC 
defined a number of sequences in the K — log P plane. Similar sequences are 
found in the Galactic Bulge (Glass & Schultheis 2003) These sequences show 
considerable internal scatter. Some of this will be due to the fact that the results 
depend on only single observations at K. However since only the Miras (in 
Wood's sequence C) are expected to have large K amplitudes (most of the stars 
are small amplitude semiregular variables) the scatter in the LMC must (except 
for the Miras) be mainly intrinsic 1 . It remains possible that such sequences 
might be narrowed by the introduction of a colour term. Even so, it is not at 
present generally possible to determine to which of Wood's sequences a field SR 
variable belongs. Thus the use of these sequences for distance determination 
seems limited for individual stars. 

7. Conclusions 

So far as can presently be determined the slope and zero-point of the Mira PL 
relation at K is invariant. The rather precise zero-point of this relation together 
with the small scatter of individual stars about the relation now makes it possible 
to derive good distance for Miras in our Galaxy and in the Local Group. The 
spectacular work of Rejkuba on Miras in Cen A indicates the potential for Mira 
studies well beyond the Local Group. Amongst important problems still to 
solve are; the evolutionary status of the hot bottom burning Miras above the 
PL relation; the nature of the Miras which apparently lie below the PL(M;, ;) 
relation in the region of the Galactic Centre; and, the nature of the Galactic 
short-period-red Miras discussed in section 3. In addition there is the problem 



In the Bulge there is addition scatter due to the depth in the line of sight. 
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of understanding the small percentage of Miras which show long term period 
changes. These have usually been attributed to shell-flashing episodes (Wood k 
Zarro 1981) but it has recently been suggested as an alternative that they are 
due to epochs of envelope relaxation (Zijlstra, Bedding k Mattei 2002). 
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